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I Abstract 

D , We examined the neutrino sector in x SU (2) p SUSY GUT with spontaneous CP violation. 

• At a glance, the discrete symmetry, which is introduced in order to solve the SUSY CP problem, 
constrains the allowed operators too strongly for the neutrino sector to be consistent with the 

CNJ ■ experimental data, i.e., the fi neutrino becomes massless as commented in the previous paper. 

^ , We showed that this issue can be solved if some operators are taken into account. And we saw 
that such operators do not play an important role in studying quark and charged lepton sectors. 

QQ , The predictions on the neutrino masses and mixings are the same as the Eq models, which are 

■ consistent with various experiments on the neutrino oscillations. 

(N 
O 
(N 

J 1 Introduction 

X 

■ Weak scale supersymmetry (SUSY) is one of the most promising candidates for physics beyond 
the Standard Model (SM) [l]. The minimal SUSY extension of the Standard Model (MSSM) has 
several attractive features. It, for instance, provides a solution to the gauge hierarchy problem and 
a dark matter candidate as the lightest supersymmetric particle (LSP). Moreover, supersymmetric 
grand unified theory (SUSY GUT) is strongly motivated by remarkable coincidence of three SM 
gauge coupling constants around 10^^ GeV. 

Grand unified theory can also explain the origin of hierarchical structures of masses and mixing 
angles in the SM particles by simply assuming that the lO's of SU (5) induce stronger hierarchical 
Yukawa structure than the 5's. However, the SUSY GUT scenario is generically suffering from 
various theoretical and phenomenological difficulties. For instance, 

1. SUSY flavor problem: if generic soft SUSY breaking terms are introduced, flavor changing 
neutral current (FCNC) processes exceeds the experimental constraints [2]. 
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2. SUSY CP problem: if generic soft SUSY breaking terms are introduced, CP violating observ- 
ables such as electric dipole moments (EDMs) exceeds the experimental constraints [2]. 

3. /i problem: supersymmetric Higgs mass /U must be the same order of the soft SUSY breaking 
scale though it is a SUSY parameter. 

4. The doublet-triplet splitting (DTS) problem and proton decay: there must be huge mass 
separation between doublet Higgs (weak scale) and triplet Higgs (above GUT scale) in order 
to make proton's lifetime longer than experimental limit. Moreover, recent proton decay 
constraint suggests that triplet Higgs mass should naively be larger than 10^^ GeV, but it 
generically spoils the success of the gauge coupling unification. 

5. Unrealistic GUT relation: generically, unification of quark and leptons in GUT tends to result 
in unrealistic Yukawa relations which are inconsistent with the observed masses and mixings 
of quark and leptons. 

Eq unification is quite attractive because it can naturally induce the feature that the lO's 
of SU{5) induce stronger hierarchical Yukawa structure than the 5's, which plays an important 
role in obtaining realistic Yukawa hierarchies in SU{5) GUT. Moreover, if we introduce a family 
symmetry, SU{3)FiSU{2)p), all three generation quark and leptons can be unified into a single 
(two) multiplet(s), and after breaking the family symmetry and the Eq symmetry, realistic quark 
and lepton masses and mixings can be obtained. Furthermore the third generation 10 of SU{5) 
can have different sfermion masses from the other sfermion masses, because it is from the third 
generation 27 of Eq. It is remarkable that such effective SUSY sfermion mass spectrum [3] can 
satisfy the LHC constraints with lighter stop, which is important to keep the naturalness in SUSY 
models. 

Unfortunately, such effective SUSY type sfermion mass spectrum is generically suffering from 
the new type of SUSY CP problem, in which Chromo electric dipole moment (CEDM) of up quark 
becomes too large to satisfy the experimental constraints. The 1-3 mixings of up-type squark masses 
become too large and have generically 0(1) complex phase after diagonalizing the complex Yukawa 
couplings if the stop mass is different from the other two up-type squarks. 

Recently a scenario which solves this problem as well as the old type SUSY CP problem by 
introducing spontaneous CP violation has been proposed The essential point is that the discrete 
symmetry, which is introduced in order to solve the old type SUSY CP problem, leads to real up- 
type Yukawa couplings which can solve the new type SUSY CP problem and to complex down-type 
Yukawa couplings which give 0(1) Kobayashi-Maskawa (KM) [5] phase. The model has remarkable 
features, for example, the doublet-triplet splitting is naturally realized by introducing the anomalous 
U{1)a gauge symmetry, predictive Yukawa matrices of up quarks, down quarks and charged leptons 
are obtained because of the reduced number of the 0{1) parameters, smaller Cabibbo-Kobayashi- 
Maskawa (CKM) [6] matrix element Vub ~ A^, where A ~ 0.22 is the Cabibbo mixing, is obtained. 

In this paper, we examine the neutrino sector in the above model with the spontaneous CP 
violation, since this has not been studied in detail in the previous paper. At a glance, the discrete 
symmetry, which is introduced in order to solve the SUSY CP problem, constrains the allowed 
operators too strongly for the neutrino sector to be consistent with the experimental data, i.e., the 
jj, neutrino becomes massless as commented in the previous paper. We will show that this issue on 
the neutrino sector can be solved if some operators, which are not important in studying quark and 
charged lepton sectors, are taken into account. And we clarify the reason why such operators play 
an important role only in the neutrino sector. The predictions on the neutrino masses and mixings 
are the same as the Eq models, i.e., sin0i2 ^ sin023 ^ X^'^, sin^is ~ A, Am^2/^"''23 ~ -^^i 
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which are consistent with various experiments on the neutrino oscillations. Since the operators are 
restricted by the discrete symmetry, the predictions on the 0{1) coefficients are expected. However, 
we are not able to find any simple predictions between the 0(1) coefficients because there are a lot 
of operators for the right-handed neutrino masses. 

As the result, we can obtain an attractive model in which 

1. realistic quark and lepton masses and mixings are obtained, 

2. all the CP phases in the model, the KM phase, CP phases in neutrino sector and CP phases in 
sfermion sectors, are given from one phase which is obtained by the spontaneous CP violation, 

3. as the result, SUSY CP problem is solved, 

4. the effective SUSY type sfermion mass spectrum is predicted. 

The paper proceeds as follows: in Section 2, we briefly review the basic features of Eq x SU (2) p x 
U{1)a_ SUSY GUT with spontaneous CP violation. In Section 3, we calculate neutrino masses and 
Maki-Nakagawa-Sakata (MNS) [7] matrix explicitly, after specifying massless 5 modes. We see that 
a specific type of higher dimensional operators is essential for the realistic neutrino masses. We also 
show that the above higher dimensional operators only affect the neutrino sector. The last section 
is devoted to summary and discussion. 

2 Basic features of Eq x SU{2)f x ^7(1)^ SUSY GUT with sponta- 
neous CP violation 

2.1 Eq unification and SU{2)p family symmetry 

First of all, 27 is the fundamental representation for the Eq group. In terms of Eq D 50(10) x U (1) v 
(and [50(10) D SU{h) x U{l)v]) it is decomposed as 

27 = 16i[10i + 5_3 + I5] + 10_2[5-2 + Sy + (1) 

where acutes are used to distinguish different 5(l)'s. Note that each 27 contains two 5's (and I's) 
of SU{5). This nature is essential for realizing different Yukawa structures of up-type quarks, down- 
type quarks, charged leptons and neutrinos from a single hierarchical structure of an Eq invariant 
Yukawa couplings [HI EJ [TO] . 

In order to break Eq gauge group into the SM gauge group GsM) we introduce three types of 
Higgs fields. First one is 27 representation Higgs H. We also introduce 27 representation Higgs H 
in order to satisfy the D-flatness condition of Eq gauge interaction. H, H obtain vacuum expectation 
values (VEVs) 

{H) = {1h) / 0, (2) 

and they break Eq into 50(10). Second one is an adjoint Higgs A, which belongs 78 representation. 
It gets Dimopoulos-Wilczek (DW) [IT] VEV proportional to U{1)b-l direction 

{A) = (45a) « U{1)b^l, (3) 

and it breaks 50(10) into 5C/(3)c7 x SU{2)l x SU{2)r x U{1)b~l- Third type are C and O, and 
they are also 27 and 27 representations. They acquire VEVs 

(O) = (16c;[9lc])/0, (4) 
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and they break SU{3)c x SU(2)l x SU{2)r x U{1)b~l into Gsm- The superpotential which gives 
superheavy masses for (5, 5) and (5, 5') pairs is 

W = Y^^i^jH + Y^^i^jC, (5) 

where {i = 1,2,3) are matter fields in family i. After H and C acquire VEVs, three degrees of 
freedom among 5j and 5^ get superheavy masses, while remaining three are massless at the GUT 
scale. We assume the Yukawa hier arches 



^.f I I (6) 

and ratio of two VEVs 
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up to 0{V) coefficients. These assumption can be realized proper charge assignment of anoma- 
lous U{1)a symmetry and the breaking effect of SU{2)f family symmetry, which we explain later 
respectively ^1. Under these assumption, three massless modes are mainly (5i, 5']^, 82). Note that 
third generation does not appear at the leading order. So we end up following different Yukawa 
hierarchies 

; \5 \3 \ / \6 \5.5 \5 \ / \7 \6.5 \6 \ 

Frf ~ I A& X* X' I , ~ Y^' ~ I A^ A^-^ A* I , ~ A" I A«-& A« A^-^ I ^^^A, 
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(8) 

when the MSSM Higgses Hu and Hd are included in IOh and 16c- Here n is a number, and we 
used the relation {C)/{H) ~ A°-^ 

In order to realize the Yukawa hierarchy between 0(1) top quark and the other quarks naturally, 
we employ SU{2)p family gauge symmetry. We treat first two generation fields as the doublet 
= 1,2) whereas the third generation fields ^'3 and all Higgs fields are treated as singlets 
under SU{2)p. Here the index a can be raised or lowered by the anti-symmetric symbols e"* and 
eab- We also introduce flavon fields Fa and F"" that are singlets under Eq and doublet and anti- 
doublet under SU{2)p, respectively. VEVs of Fa and -F"^ break SU{2)p, then this effect generate 
hierarchical Yukawa structure in Q. Since massless 5's do not contain the third generation at the 
leading order in Eq GUT, both left-handed and right-handed components of the top quark belongs 
SU{2)p singlet, while all other quarks and charged leptons are (partially or completely) SU{2)p 
doublets. 

SU{2)p family symmetry is also useful for solving the SUSY FCNC problem. Suppose that the 
soft SUSY breaking terms are mediated to the visible sector above the scale where Eq and SU (2) p 
symmetries are respected, such as in gravity mediation. Then the symmetry guarantees sfermion 
masses which degenerate at the leading order except for IO3 jlOj : 

/ ml 

mio = \ ml \, ml={ ml \ , (9) 
i / V m2 



ml 











ml 











m- 



^Strictly speaking, ([7]) is not true in U{1)a framework. The correct relation is (C) / {H) ~ A"'^, where c and 
h are U{1)a charges of C and H, respectively. 
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where mio and mg are squared sfermion mass matrices at the GUT scale in 10 and 5 fields of 
SU{5), respectively. For sfermions in 10, FCNC observables provide stringent constraints mainly 
on the first two generations, while constraints for 5's are stringent for all three generations. This 
is because the mixing angles of 10 fields are small, while those of 5 fields are large. Therefore, 
in this model, FCNC constraints can be evaded by raising mJ*^. Moreover, the weak scale is not 
destabilized as long as ms is around the weak scale. So this type of sfermion masses (" effective 
SUSY" or ''modified universality") [M El El IISl IS] solves the SUSY FCNC problem without 
spoiling naturalness. We also comment that LHC constraints to this type of sfermion masses is 
much weaker than to the constrained MSSM (CMSSM) type sfermion masses jl7] . 

Unfortunately, there is a tension between this modified universal sfermion masses and complex 
Yukawa couplings which are important to obtain the non-vanishing KM phase. In the basis in which 
Yukawa matrices are diagonal, sfermion mass matrices become 

/ A6 A3 \ 

Aio = vIomloVw - rhlo + ("^1 " "^o) A^ A^ , Ag = K^mlVg = ^g, (10) 

V A3 A2 1 / 

where Vio and Vg are the diagonalizing unitary matrices for 10 and 5 fields of SU{5), respectively, 
and we take the CKM like matrix as Viq. Since the mass insertion matrices (5io = Aio/mQ do 
not vanish in the limit mo oo, the CEDM of up quark from this SUSY contribution are not 
decoupled. The constraints for these parameters by mercury (neutron) become [13 El El [20] 

Im[(5,Ji3(5nj3i] < 3 X 10-^(9 X 10-^) {-^^^y, (11) 

which are much smaller than the prediction A^ ~ 10~^. This is a serious problem on the "effective 
SUSY" or "modified universality" , which we call the new type of SUSY CP problem. 



2.2 Anomalous U{1)a symmetry 

An Anomalous U{1)a symmetry [211 [22] is introduced in order to solve the doublet-triplet splitting 
problem and problem, and to provide the origin of the hierarchical Yukawa structures [23] . This is 
a gauge symmetry whose anomalies are canceled by the Green-Schwarz mechanism [2jj . The theory 
possesses the Fayet-Iliopoulos term and we assume its magnitude as ^ = AA. Here A is the cutoff 
scale of the theory and we set A = 1. Let us denote the gauge symmetries of the theory except for 
U{1)a as Ga- Consider a theory consisting of all the Ga and U{1)a invariant terms including non- 
renormalizable operators. Then it is shown in |23j that the theory has a supersymmetric vacuum 
where all the fields that are negatively charged under U{1)a get VEVs in the following wa\l*3|. 

(Z+)=0 (z+>0) 

(Z-)~A-^' (zr<0) ^'"^ 

Here Zj is Ga singlet field whose U{1)a charge is Zj. This argument can be extended to the case 
where Zi is composite operator that is made by Ga nonsinglet fields. For example, in Z~ = XX 
case, {XX) ~ A~(^+^) leads to \{X)\ = \{X)\ ~ X-i^+s)/'^^ once the D-flatness condition of Ga is 



*^There is upper limit of mo comes from charge and color breaking (CCB) effect. For detail discussion, see [12) . 

*^From now on, each superfield is denoted by an uppercase letter, whereas the corresponding lowercase letter 
indicates an associated U{1)a charge. The consistency of (|12p requires the number of positively charged fields to be 
larger than that of negatively charged fields by one. 
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taken account. Note that the above results can be obtained under a natural assumption that all the 
interactions which are allowed by the symmetry of the theory are introduced with 0{1) coefficients. 

It is important to mention that a term whose total U{1)a charge is negative does not appear at 
the U{1)a breaking vacuum. The reason is that this type of term should originally be accompanied 
by at least one positively charged field but its VEV is always vanishing according to (fT2]) (SUSY- 
zero mechanism) |23l I25| . A U{1)a symmetry and its specific SUSY vacuum are applied in several 
aspects of phenomenological model building. For example, an appropriate U{1)a charge assignment 
for the Higgs sector can ensure the DTS via the DW mechanism [23] , and the SUSY zero mechanism 
can be applied to solve the /i problem [26] . 

For the following arguments, we will briefly remind the points on the solution for the fi problem. 
Since the MSSM Higgses, and Hd, have negative U{1)a charges, the mass term is forbidden by 
the SUSY zero mechanism. With the positively charged singlet field S, the term SHuHd can be 
allowed, but {S) = in the SUSY limit. However, if SUSY is broken at the weak scale Ayi/, the 
VEV of S becomes non-vanishing and the order of the weak scale. This results in the SUSY Higgs 
mass with 0{K\y)- The h parameter also becomes the weak scale. 

2.3 Spontaneous CP violation and discrete symmetry 

As we have seen, one of the attractive features of x SU (2) p SUSY GUT is the modified universal 
sfermion mass (MUSM) (|9|), which evades FCNC constraints with keeping naturalness. However, 
such type of mass spectrum is generically suffering from the CEDM constraint of the up quark 
[20] . Fortunately, we can evade this constraint by introducing a spontaneous CP violation (SCPV) 
[27^128]. By using SCPV, we can realize real up-type Yukawa couplings together with the KM phase. 

Let us introduce the Eq x SU{2)p singlet field S {s > 0). Then we obtain the following 
superpotential made of flavon fields Fa^F"" and S: 

W = X'S^ c„A(-^+-^>(FF)"]. (13) 

n=0 

Here are real 0{1) coefficients and rif \s the number where SUSY-zero mechanism truncates the 
sum. When nj > 2, the F-flatness condition with respect to S leads to complex VEV (FF) and 
then CP symmetry can be spontaneously broken. Using SU{2)f gauge symmetry and its D-flatness 
condition, we can take only (F) is complex without loss of generality: 

- ( e.X-'f^m ) ' (^") - ( A-(/°+/)/^ ) • 

Unfortunately, the SCPV affects fi generation in U{1)a framework and leads to an unwanted 
outcome [4] 

Avgifib*] = 0{1), (15) 

if SFF exists. This rephasing invariant phase induces EDMs of quarks and leptons, and give rise to 
the SUSY CP problem again. In order to evade such unwanted relation, we introduce an additional 
discrete symmetry Zq, which forbids SFF. Then we can avoid this old- type SUSY CP problem. 
Moreover, this symmetry results in the real up- type Yukawa matrix, which solves new- type of SUSY 
CP problem, i.e., the diagonalizing unitary matrices become real and therefore the mass insertion 
matrices for up- type quark become real, if the MSSM Higgs Hu is included in 10 h- This symmetry 
constrains the model and reduces the number of 0(1) parameters. As a byproduct, the up-quark 
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Table 1: Field contents and charge assignment under EqxSU{2)p x U{1)a x Zq. 



Yukawa coupling y„ and the CKM matrix element Vub become closer to experimental values y« ~ A , 
Vub ~ than the expectation of Eq GUT y„ ~ A^, Vub ~ A^ [4j. 

It is also reported in that the discrete symmetry makes one of the neutrinos massless. This 
result seems to be unrealistic, but we will show that this issue can be solved by introducing some 
operators, which have not been considered in the previous paper. 

3 Calculation of neutrino mass and MNS matrix 

3.1 Field contents and Yukawa couplings 

First, we summarize the field content of the model and its representations under the Eq, SU(2)p, 
U{1)a and Zg symmetries. 

We introduce the following fields, which are listed in Table 1. All matter fields ^ belong to 27 
representation of Eq. We make its first two generations as a doublet (a = 1,2) and the third 
generation as a singlet 'I'3 of the SU{2)f, respectively. The SU{2)p and CP are simultaneously 
broken by the VEVs of the flavon fields Fa and F"". All the other fields are singlets under SU{2)f 

H is the field whose VEV (1^) 7^ breaks Eq into 50(10), and H is introduced to maintain the 
D-flatness condition. The VEV of A breaks 50(10) into SU{2>)c x SU{2)l x SU{2)r x U{1)b-l- 
C is the field whose VEV breaks 5C/(3)c x SU{2)l x SU{2)r x U{1)b-l into Gsm, and C is also 
introduced to maintain the corresponding D-flatness condition. Basically, the F-fiatness conditions 
of the positively charged fields, A' , C', C' etc., determine the VEVs of the negatively charged fields. 
For example, the F-flatness conditions with respect to A' make A acquire DW-type VEVs for the 
50(10) adjoint component to solve the DTS problem. The alignment between the VEVs of A and 
C, C are realized and at the same time, the pseudo Nambu-Goldstone modes become heavy [29]. 
In Table 1, U{1)a charges are assigned so that the DTS and appropriate Yukawa hierarchies are 
realized. Also, Zq charges are determined so that the SUSY CP problem is evaded. 

5, S' and S" are introduced to realize the SCPV and /i generation. The relevant superpotential 

To cancel the Witten's anomaly, odd number of additional doublets of SU (2) f are required. 
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are 



Ws = \'S + X'+^^^SZl + X'+^''SH^, (16) 
Ws' = A^'+^2y^2 + A"'+^+"5'C'C + A^'+2/+2/5'(FF)2, (17) 
Ws" = X'"+'^S"Z2 + X'"-^~'^'S"CC. (18) 

First, we obtain the VEV {Z2) from the F-flatness condition of Ws with respect to S. Since the 
one of the three solutions is real, we assume (Z2) is real. Second, the F-flatness condition of Ws" 
determines the real VEV (CCj i. Third, the F-flatness condition of Ws' leads the complex VEV 
(FF), which breaks the SU{2)p and CP spontaneously. Note that the operator gives the SM 
Higgs mass after developing the VEV of H if the SM Higgses are included in IOh- Then, after 
breaking SUSY, Ws generates fi and b terms. Since does not couple to S" and S", so fi and b 
are real. 

Let us examine the mass matrices of quarks and leptons in this model. Under the charge 
assignment of Table 1, the following interactions between matter and Higgs fields are allowed: 

/ d^''{A,Zs,HH)^a \ 

: d^''{A,Z^,HH)^a cX^^^^+^^F'^maF^^b bX^t'^+^^+f F"^ ^ \>^'^H, (19) 
V hX^-+'^'^+f^^F''^a aX^^-'^^i^z ) 

/ fX^^^+f+^F^^aF^^h g'X^'-+'^-'+fF''ma^z \ 

: f'X^^-+f+fF^^aF''^b X^C. (20) 

\ g'X'^-+^--'+f'^3F''^a / 

Here we explicitly write in all 0{1) coefficients, o, b, c, d, f and g' . Note that all the 0(1) coefficients 
are assumed to be real because of the original CP symmetry. 

In (llOp . the structure of (1,2) and (2,1) elements is a little bit complicate, so we should see 
this point carefully. ^""{A, Z^, HH)^a consist of several types of terms. Cleally, some additional 
field(s) for (1,2) and (2,1) elements are needed because ^"■^a = e°-^^a^b = 0, where e^^ = -e^^ = 
l,e^^ = e^^ = 0. The negatively U{1)a charged fields which can have non-vanishing VEVs can be 
a candidate for the additional fields. There are two possibilities. One is using the adjoint Higgs A: 
^"AZs^a, '^"A'^^a, etc. Since the VEV of A breaks 50(10) into 5C/(3)c x SU{2)l x SUi2)R x 
U{1)b-l, the contributions of these terms are different for the different components of under 
SU{3)c X SU{2)l X SU{2)r x U{1)b-l- Because the VEV of ^ is proportional to the B - L charge, 
there is no contribution for the field which has vanishing B — L charge. The other possibility, which 
has not been considered in the previous paper, is using HH: ^"^HH^aj etc. Since the VEV of 
H breaks Eq into 50(10), it is useful to write down the contributions with the representations 
of 50(10). For example, e'^^lQ^^l^,^ and e"*16*„10^,. Note that the terms e^^'lOvi-^lO^, and 
e"^16>ii^l6.lij are vanishing and give no contribution. Here we use a conventions for the component 
fields of 50(10) of a field X as 

X = 27x = 16x + lOx + Ix, (21) 
16x = Qx + Ux+Ex+Dx+Lx+Nx, (22) 
IQx = H^ + {Hu)x + H^ + {Hd)x. (23) 



*^We must forbid the term S< TFCC, which lead the complex VEV (CC). If we can choose the basis in which 
(C) is real and {0) has an opposite phase to {F), the up-type Yukawa can keep real. However, it changes the phase 
of the down-type Yukawa into removable, and the KM phase cannot be realized. 
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Then, we can write down the concrete components of the interactions ^'^(^, ^3, HH^^a 

where d^, dq, di, h are the 0{1) coefficients that are different from each other generically. As we 
will discuss later, the h term is very important for neutrino masses. 



3.2 Massless modes 



Once we write down the Yukawa couplings, we can calculate the massless combinations of 5j and 
5-. We mainly follow the procedure given in [30]. First, we fix 0(1) coefficient of two VEVs: 



X^^ 



xX 



0.5 



(25) 



where x is a real 0(1) coefficient. Then, after the Higgs and fiavon fields acquire VEVs, (|19|) and 
(|20p induce following mass matrix for 5,, and (5^,5j): 








ad^X^ 
cX^ 
6A2 





6A2 

a 



fe'PX^-^ 
fe'PX^-^ 
ge'PX^-^ 




{Ml M2 



(26) 



Here, we reparametrize 0(1) coefficients as / = xf and g = xg' . Each power of A is determined 
by the corresponding U{1)a charge. It is important to note that a = 1 for the colored Higgs 
components (i?*^, H^) of 5 and 5', and a = for the doublet Higgs components (//„, H^) of those, 
since (1, 2) and (2, 1) elements of (I26p originate from the B — L conserving VEV of A. 

In order to find the massless combinations, let us diagonalize the 3x6 matrix ()26p as follows: 



V\Mi Ms) 
Here ^ is a 3 x 3 unitary matrix and 

U 



^10 
^16 



^10 
^16 



(a4'"s 0). 



^10 



^1% 



(27) 



(28) 



is a 6 X 6 unitary matrix, and they rotate 5, and (5^,5j) (i = 1,2,3) into their mass eigenstates, 
respectively. Our task is finding matrices and C/{g which are related by 



MiC/^o + MsC/^e = 0. 
According to the calculation in Appendix A, they become 



f^io 



aadsibg-af) \ 2.5 ^ip 



ip 



ac—W 

\ (9 \b bg-af \ \3.5„ 



1 

aad; 



5 X 



ac—b'^ ' 
bads \3 
ac—b'' 



0(A5-5) 
0(A4-5) 
0(A6-5) 



1 

0(A6) 

bg-af ads ^3 
ac—h'^ g 






""'i a \2.5, 
'TiZ^g^ ^ 



-ip 




(29) 



(30) 



(31) 
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The calculation is taken at the leading order, but (C/iq) 13,23,33 ^^'^ (^^16)21 ^^'^ calculated at the 
next-leading order because the corresponding leading terms are cancelled. Therefore, the massless 
modes of 5 can be written as 

5° = (C/°t) . .5^. + (f/Ot) . .5^. ^1 5/ + . . . I . (32) 

In order to calculate neutrino mass, we should also specify the MSSM Higgs doublets and Hd- 
They are the combinations which remain massless at the GUT scale. According to the Appendix 
A in they are 

Hu C 5h (33) 

HdC5'jj + f3He-'fX'>-''5c. (34) 

Here is a real 0(1) coefficient and the phase e~^P comes from VEVs of the F and As 
mentioned in [3] , the second term in p4p is important to ensure the KM phase and make all charged 
leptons massive. 

3.3 Neutrino mass 

Now we calculate neutrino masses in our model. In this subsection we explicitly write the cutoff 
scale A. First, we compute the Dirac neutrino mass matrix. Since we have 6 right-handed neutrino, 
Ij and 1[, from the interactions lOjl^lO// and 16jl6jl0//, 3x6 Dirac neutrino Yukawa matrix 
Yjy^ is given by 

= {Y^^v = {U?^{Y,^f U'lliYfn (35) 



where 



It is also important that a = in the mixing matrices (|30|) (j31|) . Using these facts, we can express 
the Dirac mass as follows: 



hX^ \ / 

= { -hX^ cA^ 6A2 , Yf=i yiX^ cA^ bX^ | . (36) 
6A2 a / \ 




h^^X^-^e'P -fX^-^e'P -gX^-^e'P 
5 



ac—TP' 







Y^^= \ -/lA^ 



-fdzA^ 

I • (37) 

IdiX^ ic-'{)X^ ^A2 



Note that the second generation would be massless without h. The contribution from the term 
^'^HH^aH is essential for obtaining the realistic neutrino sector. 

Second, we introduce the 6x6 Major ana right-handed neutrino mass matrix. Since each 27 has 
two I's, the Majorana mass has following three types of contribution: 



+ + ^ia,(lc)(lc)- (38) 



®If we choose Z2 = —1 (or 22 = —2), then the term C' Z2C (or C"Z|C) changes the phase in (|34p into arbitrary. 
For detailed discussion of the Higgs sector, see the Appendix A in [4]. 
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Then the Majorana mass term is 6 x 6 matrix 
Here, we parametrize with 0{1) coefficients as 



Y,l Y 



11 



Y, 





-13 


iVl2A' 


-12 


NrsX- 


-10 


NuX~ 


12.5 


iVi5A" 


11.5 


NieX- 


9.5 


N12X 


-12 


iV22A- 


-11 


N23X 


-9 


N2iX~ 


11.5 


iV25A^ 


10.5 


N26X- 


8.5 


N13X 


-10 


A^23A 


-9 


N33X 


-7 


N3aX- 


-9.5 


iVssA- 


-8.5 


N36X- 


6.5 




V2.b 


iV24A- 


11.5 


N3aX- 


9.5 


N^aX^ 


-12 


iV45A- 


-11 


N^eX 


-9 


iVisA- 


11.5 


iV25A- 


10.5 


N35X- 


8.5 


Ni5X' 


-11 


iV55A- 


-10 


N^eX 


-8 


iVieA- 


-9.5 


N2eX'- 


-8.5 


iVaeA" 


6.5 


N^eX 


-9 


iV56A 


-8 


NqgX 


-6 



J 

where Nij are 0{1) complex parameters. 

Finally we can calculate the neutrino mass using the seesaw mechanism. The result is 

yiiA"^ yi2A"^-^ yisA^^ 



ITT \2 

M„ = Y„ -^-^ , = Y„^Y^}Y,T 



A 



yi2A-i-5 
yisA"^ 



y22A 2 y23A 
y23A"^-^ yssA"^ 



where yij are complicate combinations of C(l) coefficient: 



yii = Nii{h 



bg - af 2 



e^'P + iV22/"e 



ac 



62 



2^2ip 



ac 



62 



+ 2N23fge^'P - 2N25f{-^di)e'P - 2N35g{-ldi)e'P , 



yi2 



bg - af 

— r 



ac 



N22fe'P - N23ge'P + N- 



251 



yi3 



h 



bg - af 



ac — 62 
fe'P 



,3 , ,ac — 

2 a 



62 bq — af b , bq — af 
— + , -) + Nie^ 



g 



N2A{^di) + N25{ 

2 a 



b a 



-ge 

3 



,T /3,, ,^ .ac — b"^ bg — afb. 

N34i-di) + A^35 + , - + N36 

2 a b a 



di 



y23 



?/33 



,x /3 , , ,^ ,ac — 

A^45 T^dz + ^^55 

2 a 



N24i-dl) + N25{ 



bg 


-an 




g 


bg 


-an 




g _ 


bg 


-af] 



6 a 
y22 = fi^ N22, 

+ bg-af b^ ^ ^r j9-af 
a b a g 



N,,{\dif + NU- 



a b a 2 q 



g 



-62 bq - af b. 

+ . -)• 

a b a 



(39) 



(40) 



(41) 



(42) 
(43) 



(44) 
(45) 
(46) 



(47) 



These masses give Amf2/Am|3 ~ A2 = 4.8 x 10 2^ which is consistent with the experimental 



facts Amf2/A 



m 



23 



3.1 X 10 |31j up to the combination of 0(1) coefficients. 
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3.4 MNS matrix 

We can also calculate the MNS matrix in our model. First, charged lepton Yukawa matrix is given 
m as 

Ye=\ ^Hfe'^X"-'' PHg'e'^X^-'' . (48) 

\^ -fd,A5 (-^ + M^A)A4 fa/A2 ) 

Then we can calculate MNS matrix. Following the diagonal procedure given in [32] , it is computed, 
at the leading order, 

/ 1 z;i2A0-5 visX \ 

Vmns = -^^i*2A°-' 1 ^^23A°-^ , (49) 

v ivt,v*,,-vi,)x -v*,,x'-' 1 ; 

where vi2 , fia , and V23 are written as 

^12 = i rbTT? 2 — ' (^0) 

Ph 1 - 2/22^33 - 

yi3 * 2/12^33 - ^132/23 .^-.x 

i;i3 = V23 —, (51) 

y33 y222/33 - ^23 

I^Hgg' y23 

V23 = 7 7 • 52 

bg - af 1/33 

This result can be compared with experimental values. The model says tan0i2 = Vvil^w ~ 
^0.5 _ Q_47^ tan 023 = V23/V33 ~ A*^'^ = 0.47 up to the combination of 0(1) coefficients. These are 
consistent with experiment tan0i2 ~ 0.68, tan 023 > 0.747 [3l]. It is also interesting that recent T2K 
result [33] 0.087 < sin0i3 < 0.275 (assuming bcp = 0) supports our result sin 6*13 = 1^13! ~ A = 0.22, 
up to the combination of 0{1) coefficients. 



3.5 The role of the ^!HH^!H terms 

As we have seen, the terms ^°^HH^ aH played the important role for the non-vanishing second 
generation neutrino mass. Here we show that this type of higher-dimensional operators do not 
affect other structures of this model, for example, the quark masses and charged lepton masses. 

First, we classify 50(10) invariant Yukawa interaction derived from Eq invariant Yukawa term 
^!i^>jH (or ^.i^jC): 

27i27^27H,c = 16il6jlOH,c + 10^10^1^,0 + 16il0jl6H,c + 10iljlOH,c- (53) 

Note that these terms have common 0(1) coefficients, and all terms are symmetric between i and 
j, although at 50(10) level, the coefficients of the terms IBilOjlBi^^c and 10iljlOH,c do not have 
to be taken as symmetric. 

This property changes if we include higher-dimensional operators like ^{HH^jH or ^iHH'^jC. 
Each terms in the right-hand-side of (j53p receives extra contribution, which generically differs from 
each other. Therefore, they cannot have common 0(1) coefficients, and the number of parameters 
increases. 

Fortunately, in our model, only (1,2) and (2,1) elements of have room of U{1)a charge for 
including HH, and in other Yukawa elements, SUSY-zero forbids such higher-dimensional operators. 
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Since these elements are anti-symmetric due to SU{2)f, the extra contributions 16j 16^10// and 
lOjlOjl/f vanish. Therefore, the higher-dimensional operator including HH only affects the Dirac 
neutrino mass through lOjljlO/f, while the quark sector and charged lepton sector is unchanged. 

According to the previous paper the model has characteristic predictions Vub ~ and 
VcbUb = He at the GUT scale. We will show in Appendix B that these predictions are not spoiled 
even if the model generically contains such higher-dimensional operators. Sufficient conditions for 
these predictions are {A) oc Qb-l and that the (2,2), (2,3), (3,2), and (3,3) components of 
are vanishing. Such conditions are satisfied in more general Eq GUTs with family symmetry. The 
former condition should be satisfied to solve the doublet-triplet splitting problem, and the latter 
can be satisfied by the SUSY zero mechanism if the anomalous U{1)a charge of C is smaller than 
that of H. The discrete symmetry for solving the SUSY CP problem is not necessary. Therefore, 
these two predictions are rather general ones. The first prediction Vub ~ A^ has been tested by 
B-factory experiments, and the second prediction VcbVb = Vc can be tested in future experiments by 
measuring the tan/3 = {Hu)/{Hd). 

4 Summary and discussion 

In this paper we conducted a detail calculation of the neutrino sector in Eq x SU{2)f SUSY GUT 
with spontaneous CP violation. Originally, Eq GUT can explain the realistic neutrino sector as 
well as quark and charged lepton sector. However, the previous paper reported the problem about 
the neutrino mass because of the new discrete symmetry, which is introduced in order to solve the 
SUSY CP problem. We computed the neutrino masses and the MNS matrix explicitly and found 
that the term like ^iHH^jH, which have not been considered in the previous paper, is important 
to make the second generation massive. These higher dimensional operators only contribute in the 
neutrino sector, so the structures of quarks and charged leptons, which calculated in the previous 
paper, are kept valid. 

Our result reproduces experimental value well. Am^g/^^'-isi taii0i2, tan023 and sin^ia are 
all consistent with experimental values up to the combination of 0{1) coefficients. Combining the 
previous paper and our result, the concrete model reproducing realistic masses and mixings for all 
the SM fermions is constructed. Note that even if the model has the "modified universal sfermion 
mass spectrum", we have no SUSY CP problem. 

In this model, because of the strong constraint of Yukawa structure, there are characteristic 
predictions Vub ~ A^ and VcbUb = Vc at the GUT scale [1] . These predictions are valid even if we 
generically include the higher-dimensional operators ^iHH^jH and ^iHH^jC. Unfortunately, we 
could get no characteristic prediction in the neutrino sector because of huge number of parameters 
in the Majorana neutrino mass terms. 
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A Diagonalization procedure of the superheavy Yukawa matrix 



In this appendix, we derive the expression for the two matrices [/{g (j30p and [/{g (j3ip which describe 
the massless 5? as the combination of 5j and 5^ as All calculation is performed at the leading 
order. 

First, we show the strategy of finding U which diagonalize the superheavy Yukawa matrix as 
()27p . In the following, we assume that Mi has rank 3, and therefore we have a inverse matrix M^^. 
Then we calculate M^^M2 and parametrize it as 



2.5 



Mr ^ A/2 



c 



C2lAl-5 C22A0-5 C23A"^-5 
,A3-5 C32A2-5 C33A0-5 



(54) 



C31/ 



Here Cij are the combination of 0{1) coefficients. After expressing C/J^q and ?7j^g in terms of Cjj, we 
will substitute concrete parameters in them. Then we rotate M^^^Mi M2) = (13x3 C) as 



(13X3 C)U' = {C 03x3), 

where ^7' is a 6 x 6 unitary matrix 

■ Ui Us 

and C" is a 3 X 3 matrix. If we find unitary matrices V and V' which diagonalize MiC' as 



U' = 



then we can express the matrix U as 
U = 



U^ r/0 

^10 "^lO 



^16 



•^16 



u' 



V 






T 



UiV U3T 
U2V' UaT 



(55) 
(56) 

(57) 
(58) 



Here T is a unitary matrix which corresponds to the degree of freedom of the rotation for the three 
O's with keeping (p7|) . 

Since we are only interested in C/^q and U^q in this paper, our task is to find the matrix U' . 
After tedious calculation, we can find C and U' as 



:i3x3 C)U' 




(C22 





1 

^)(^32 




u' 



( ^A2-5 
' Ci3 




Cl3 

Cl2 \2 



C13 
1 





-(C22 



.£33A3 

Cl3 



C13 
1 



0(C32-)A^ 



(C21 - 

(C22 - 

- — (C22 
Ci3 V 



C11C23 ) \T75 

C13 

C12C23 \ \0.5 
C12C23 \ \2.5 



(C31 
(C32 



■tt(-32 



CllC33 \ x3.5 
C13 

C12C33 \ 1, 2.5 

_ Ci2C33\^4.5 
Cl3 



- ^)A2-5 


1 







-(C2i - ^)Ai-5 

-(C31 - ^)A3-5 


£23A 

C33A3 

CI 3 


-(C22 
-(C32 


C12C23 j -yCS 

Cl2'^C33 ^ \2.5 
C13 M 


1 










-(c21-^)(c22-^)A2 
_ciiA3 

Cl3 


ff(c22-^)Al-^ 

_^A2-5 

Cl3 




1 

_£i2A2 

C13 



(59) 



(60) 
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We choose a basis which (C/iq) 12, (f^i6)i3 ™d (?7ig)22 are zero by the rotation 

/I \ 

r= 1 ^(c22-^)Ai-5 . (61) 

\o -^(C22-^)A^-^ 1 y 

After that, we can get [/{g and Uiq in terms of Cij as 

ft(c21 - ^-^)X'-' 1 ft(c22 - ^)A^-5 \ 

U?o = I -(C2i - ^)Ai-^ fff A -(C22 - ^^)A0-5 , (62) 

-(c31-^)A3-5 -^A3 (c32-^^)A2-5 / 

1 \ 

<=| -(C2i - ^)(c22 - ^)A2 1 . (63) 

_ai;^3 _J__\2.5 _£i2;\2 / 

Cl3 Ci3 C13 / 

We can check these matrices satisfy the relation ([29l) at the leading order. Finally, we obtain (pOj) 
and (j3T]) by substituting the concrete expressions in Cij: 

bg-af ac-b'^f ac-b'^g 
Cll = 1 , C12 = jj ' '^13 = 32 (6^) 

f g 

C21 = 0, C22 = —r^ C23 = —3- (65) 
ads ads 

g bf bg 

C31 = C32 = C33 = — (66) 

a aaa^ aaa^ 

B The condition of the predictions 

In section 3, we obtained the realistic neutrino sector by introducing the higher-dimensional operator 
^'^HH^ aH , and saw that such operators cannot be written in other elements of and because 
of the SUSY zero mechanism. In this appendix, we will clarify the conditions of the characteristic 
predictions V^h ~ A^ and VcbUh = Vc, and show that the conditions can be satisfied even if we 
generically include the higher-dimensional operators ^iHH^jH and '^iHH^jC. 

In order to see the essence of the predictions, let us choose the basis where the up-type Yukawa 
matrix Y^ is diagonal. In this basis, Vub is simply written only by the down-type Yukawa matrix 
as Vub ~ (^d)i3/(^d)33- The relation VcbVb = Vc is also simplified as {Yd)2?. = {Yu)i2- Therefore, we 
focus on the quantities (Id) 13 and {Yd)2?,- 

Let us determine the third generation of massless 5 combination 83, which is mainly 82, in this 
basis. Since we can choose the basis where 83 does not contain 5i and 5'^, we parametrize 

As we saw in Section 3.2, the massless modes are determined by diagonalizing the 3x6 mass matrix 
(Ml M2). Both Yu and Mi is obtained from Y^ , but here we consider general case where these two 
matrices cannot be diagonalized simultaneously. In this case. Mi is the form 

1 A A3 \ / a6 A5 A3 
Ml - I A 1 A2 A5 A^ A2 I I A 1 A^ I - I A5 A4 A^ I . (68) 
A3 A2 1 / \ A3 a2 1 
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On the other hand, M2 is obtained from Y'^ . Again, we consider general case 



Mo 





^5.5 

^4.5 ^2.5 
^2.5 ;^0.5 



(69) 



Then we can obtain the massless combination of 63 by solving 

/ \ 

a 

/? 

1 

V 7 y 



(Ml M2) 



(70) 



The solution becomes generically 



a ~ 



2,5 



7 



(71) 



Let us evaluate (Yd) 13 and (1^)23 to obtain the conditions for Vub ~ A^ and VcbVb = Uc- The 
down-type Yukawa matrix consists of two types of interaction: one is the interaction between 
10 and 5 via Y^^^, the other is between 10 and 5' via 

;y0.5yC_ jjere we denote the matrix Y 
appearing in the down- type quark Yukawa matrix as Yf^ . Y^ is written in this basis as 



1 ~ 



1 
A 



A3 A 




A5 A^ A2 
A^ A3 A 
A2 A A-i 




A5 A4 A2 
A^ A3 A 
A2 A A-i 



(72) 



For example, the leading contributions for (1^)13 are {Y^^Ji2, 7(^iLn)i3; aX^-^(Y'^)i2, and 
(3X^'^{Y^)i^. They are of order A^, and therefore Vub becomes 0{X^). In order to obtain Vub ~ A^, 
all the leading contributions must vanish. The first two contributions vanishes if (A) is proportional 
to Qb-l, which plays an important role in solving the doublet-triplet splitting problem. This is 
because Y^^^ is the same as Yu and therefore off-diagonal elements of Y^^^ are vanishing in the 
basis in which Yu is diagonal. Another important point is that the U{1)a charge of C is smaller 
than that of H by one. This nature forbids the (2,2), (2,3), (3,2) and (3,3) elements (before rotating 
the basis) of Y^ by the SUSY -zero mechanisnO- Therefore, after rotating, Y^ becomes 



Y 



c 




A5 A^ A2 
A4 a5 A3 
A2 A3 A5 



(73) 



This changes the determination of massless 53 composition ()7ip into 



a ~ A^-S, f3 ^ A*'^ 7 ~ A 



4,5 



(74) 



so the terms containing a and /3 become smaller. The above two points are the essential points for 
Vub = at the leading order. 



*^In our model, the (1,1) entry of is also forbidden by the discrete symmetry. 
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Next, we see (1^)23- The main contributions are from 7(^ckLn)23, (^iLn)22, aA°'^(^'^)22 and 
^A°'^(y'^)23, but the first contribution is vanishing if (^4) is proportional to The second 

term is nothing but (1^)22) so the relation VcbUb = Uc is realized when a and /3 become smaller due 
to the SUSY-zero mechanism as noted above. 

How does above discussion changed if we generically include the higher-dimensional operators 
'^iHH^ jH and '^iHH'^ jHl First, Yu and ^^^^n ^"^^ derived from same expression 16jl6jl0jy 
at 50(10) level, therefore simultaneous diagonalizability of and Y^^^^ does not affected by the 
higher-dimensional operators. Second, the form of Y^ with the SUSY-zero mechanism (j73p does 
not changed because adding HH (h + h < 0) cannot revive the terms which was forbidden by the 
SUSY-zero mechanism. Therefore, the characteristic predictions Vub ~ A"^ and VcbVb = Uc do not 
affected whether or not we include such higher-dimensional operators. 

What becomes obvious in this appendix is that the conditions for Vub ~ and VcbVb = Uc can 
be satisfied in the more general Eq GUTs with family symmetry. The discrete symmetry for solving 
the SUSY CP problem is not necessary. One prediction Vub ~ A^ has been already tested by the 
B-factory experiments. The other prediction VcbVb = Uc at the GUT scale means that comparably 
small tan/3 = {Hu) / {Hd) ^ which can be tested in future experiments. 
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